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ABSTRACT 

When building a cylindrical algebraic decomposition (CAD) 
savings can be made in the presence of an equational con¬ 
straint (EC): an equation logically implied by a formula. 

The present paper is concerned with how to use multiple 
ECs, propagating those in the input throughout the pro¬ 
jection set. We improve on the approach of McCallum in 
ISSAC 2001 by using the reduced projection theory to make 
savings in the lifting phase (both to the polynomials we lift 
with and the cells lifted over). We demonstrate the benefits 
with worked examples and a complexity analysis. 

Categories and Subject Descriptors 

1.1.2 [Symbolic and Algebraic Manipulation]: Algo¬ 
rithms— Algebraic algorithms, Analysis of algorithms 

General Terms 

Algorithms, Experimentation, Theory 


Keywords 

cylindrical algebraic decomposition, equational constraint 


1. INTRODUCTION 

A cylindrical algebraic decomposition (CAD) splits R" into 
cells arranged cyhndrically, meaning the projections of any 
pair are either equal or disjoint, and such that each can be 
described with a finite sequence of polynomial constraints. 

Introduced by Collins for quantifier elimination in real 
closed fields, applications of CAD include: derivation of op¬ 
timal numerical schemes 18 , parametric optimisation [19| , 
epidemic modelling , theorem proving , reasoning with 
multi-valued functions [^, and much more. 

CAD has complexity doubly exponential in the number of 
variables [^. For some applications there exist algorithms 
with better complexity (see [^), but CAD implementations 
remain the best general purpose approach for many. This 
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may be due to the many extensions and optimisations of 
CAD since Collins including: partial CAD (to lift only when 
necessary for quantifier elimination); symbolic-numeric lift¬ 
ing schemes |29[|22| ; local projection approaches [^[^; and 
decompositions via complex space [11[ [^. Collins original 
algorithm is described in while a more detailed summary 
of recent developments can be found, for example, in [^. 

1.1 CAD computation and terminology 

We describe the computation scheme and terminology that 
most CAD algorithms share. We assume a set of input poly¬ 
nomials (possibly derived from formulae) in ordered vari¬ 
ables X = xi -< ... < x„. The main variable of a polynomial 
(mvar) is the greatest variable present under the ordering. 

The first phase of CAD, projection, applies projection op¬ 
erators repeatedly, each time producing another set of poly¬ 
nomials in one fewer variables. Together these contain the 
projection polynomials used in the second phase, lifting, to 
build the CAD incrementally. First R is decomposed into 
cells which are points and intervals according to the real 
roots of polynomials univariate in xi. Then R^ is decom¬ 
posed by repeating the process over each cell with the bi¬ 
variate polynomials in (xi,X 2 ) evaluated at a sample point. 

This produces sections (where a polynomial vanishes) and 
sectors (the regions between) which together form the stack 
over the cell. Taking the union of these stacks gives the CAD 
of R^ and this is repeated until a CAD of R" is produced. 

At each stage cells are represented by (at least) a sample 
point and an index. The latter is a list of integers, with 
the fcth describing variable Xk according to the ordered real 
roots of the projection polynomials in ( 2 : 1 , ... ,Xk). If the 
integer is 2i the cell is over the ith root (counting from low 
to high) and if 2i -|- 1 over the interval between the ith and 
(i -I- l)th (or the unbounded intervals at either end). 

The projection operator is chosen so polynomials are de- 
lineable in a cell: the portion of their zero set in the cell 
consists of disjoint sections. A set of polynomials are de- 
lineable if each is individually, and the sections of different 
polynomials are identical or disjoint. If all projection poly¬ 
nomials are delineable then the input polynomials must be 
sign-invariant have constant sign in each cell of the CAD. 

1.2 Equational constraints 

Most applications of CAD require truth-invariance for log¬ 
ical formulae, meaning each formula has constant boolean 
truth value on each cell. Sign-invariance for the polynomi¬ 
als in a formula gives truth-invariance, but we can obtain 
the latter more efficiently by using equational constraints. 



Definition 1. We use QFF to denote a quantifier free 
Tarski formula: Boolean combinations (A,y,~^) of state¬ 
ments about the signs (— 0, > 0, < 0) of integral polynomials. 

An equational constraint (EC) is a polynomial equation 
logically implied by a QFF. If an atom of the formula it is 
said to be explicit and is otherwise implicit. 

Collins first suggested that the projection phase of CAD 
could be simplified in the presence of an EC [^. He noted 
that a CAD sign-invariant for the defining polynomial of an 
EC, and sign-invariant for any others only on sections of that 
polynomial, would be sufficient. An intuitive approach to 
produce this is to consider resultants of the EC polynomial 
with the other polynomials, in place of them. This approach 
was first formalised and verified in |24| . 

A recent complexity analysis showed that using an EC 
in this way reduces the double exponent in the complexity 
bound for CAD by 1. A natural question is whether this can 
be repeated in the presence of multiple ECs. An algorithm 
for CAD in the presence of two ECs was detailed in [25| . 
The main idea was to observe that the resultant of the poly¬ 
nomials defining two ECs is itself an EC, and so the same 
ideas could be applied for the second projection as for the 
first. However, this approach was complicated as the key 
result verifying could not be applied recursively. 

1.3 Contribution and plan 

This paper discusses how we can extend the theory of ECs 
to produce CADs more efficiently. In Section [2.1| we revise 
key components of the theory for reduced projection in the 
presence of an EC from [24[ |25| . Then in Section |2.2| we 
explain how it can also give reductions in the lifting phase, 
allowing us to propose and verify a new algorithm in Section 
[3] for making use of multiple ECs. This breaks with the tra¬ 
dition of producing CADs sign-invariant for EC polynomials, 
instead guaranteeing only invariance for the truth of their 
conjunction. We demonstrate our contributions in Sections 

and 1^ with a worked example and complexity analysis. 

All experiments in Maple were conducted using Maple 
18. All code and data created for this paper is openly avail¬ 
able from http : //dx. doi . org/10.15125/BATH-00071 

2. CAD WITH MULTIPLE EQUATIONAL 
CONSTRAINTS 

2.1 Key theory from | [23|[24|[2^ 

We recall some of the key theory behind McCallum’s oper¬ 
ators. Let cont, prim, disc, coeff and Idcf denote the content, 
primitive part, discriminant, coefficients and leading coeffi¬ 
cient of polynomials respectively (in each case taken with 
respect to a given mvar). Let res denote the resultant of a 
pair of polynomials. When applied to a set of polynomials 
we interpret these as producing sets of polynomials, e.g. 

res(A) = {res(/i, fj) \ fi £ A, fj £ A, fj yf fi} . 

Recall that a set A C Z[x] is an irreducible basis if the el¬ 
ements of A are of positive degree in the mvar, irreducible 
and pairwise relatively prime. Throughout this section sup¬ 
pose B is an irreducible basis for a set of polynomials, that 
every element of B has mvar Xn and that F C B. Define 

P(B) := coeff(B) U disc(i3) U res(i3), (1) 

Pf(B) := P{F) U {res(/, 5 ) \f£F,g£B\F}, (2) 


123,24 25 


Pp{B) ~ Pf{B) U res(B \ F), 


( 3 ) 


as the projection operators introduced respectively in [23[ 
2^ |25| . In the general case with A a set of polynomials 
and E A we proceed with projection by: letting B and 
F be irreducible basis of the primitive parts of A and E 
respectively; applying the operators as defined above; and 
then taking the union of the output with cont (A). 

The theorems in this section validate the use of these oper¬ 
ators for CAD. They use the condition of order-invariance, 
meaning each polynomial has constant order of vanishing 
within each cell, which of course implies sign-invariance. We 
say that a polynomial with mvar Xk is nullified over a cell 
in if it vanishes identically throughout. 


Theorem 1 ( [23| ). Let S be a connected submanifold 

o/R”“^ in which each element of P{B) is order-invariant. 

Then on S, each element of B is either nullified or analytic 
delineable (a variant on delineability, see \2SI( ). Further, the 
sections of B not nullified are pairwise disjoint, and each 
element of such B is order-invariant on such sections. 


Suppose we apply P repeatedly to generate projection poly¬ 
nomials. Repeated use of Theorem [^concludes that a CAD 
produced by lifting with respect to these projection polyno¬ 
mial is order-invariant so long as no projection polynomial 
with mvar Xk is nullified over a cell in the CAD of (a 
condition known as well-orientedness which can be checked 
during lifting). If this condition is not satisfied then P can¬ 
not be used (and we should restart the CAD construction 
using a different projection operator, such as Hong’s 20 ). 


Theorem 2 ( [24| ). Let f and g be integral polynomials 

with mvar Xn, r(a;i,..., a;„_i) be their resultant, and sup¬ 
pose r yf 0. Let S be a connected subset ofW^~^ on which f 
is delineable and r order-invariant. 

Then g is sign-invariant in every section of f over S. 


Suppose A was derived from a formula with EC defined by 
E = {/}, and that we apply Pe{A) once and then P repeat¬ 
edly to generate a set of projection polynomials. Assuming 
the input is well-oriented, we can use Theoremj^to conclude 
the CAD of R"“^ order invariant for Pe{A). The CAD of 
R" is then sign-invariant for E using Theorem and sign- 
invariant for A in the sections of E using Theorem]^ Hence 
the CAD is truth-invariant for the formula. 


What if there are multiple ECs? We could designate one 
for special use and treat the rest as any other constraint 
(heuristics can help with the choice [^). But this does not 
gain any more advantage than one EC gives. However, we 
cannot simply add multiple polynomials into E at the top 
level as this would result in a CAD truth-invariant for the 
disjunction of the ECs, not the conjunction. 

Suppose we have a formula with a second EC. If this has 
a lower mvar then we may consider applying the reduced 
projection operator again at this lower level. In fact, even 
if the second EC is also in the mvar of the system we can 
propagate it to the lower level by noting that the resultant 
of the two ECs is itself an EC in one fewer variable. 

So we consider applying first the operator Pe{A) where 
E defines the first EC and then Pe'{A') where A' = Pe{A) 
and E' C A! contains the EC in one variable fewer. Unfor¬ 
tunately, Theorem does not validate this approach. While 
it could be applied once for the CAD of R"“^ it cannot then 









Figure 1: The polynomials from Example [l] 


validate the CAD of R" because the first application of the 
theorem provided sign-invariance while the second requires 
the stronger condition of order invariance. Note however, 
that this approach is acceptable if n = 3 (since in two vari¬ 
ables the conditions are equivalent for squarefree bases). 

Example 1. The following are graphed in Figure^ with 
g the sphere, /i the upper surface and f 2 the lower: 

p I 2 i p 2 I 2 I 2 I 2 -1 

fi = X + y + z, f 2 = X - y + z, g = x+y+z-l. 

We consider the formula (j) = f\ = 0Af2 = 0Ag>0. The 
surfaces fi and f 2 only meet on the plane y — 0 and this 
projection is on the right of Figure^. From this it is clear 
the solution requires |a:| > \/2/2 and z = —x. 

How could this be ascertained using CAD? With variable 
ordering z y x a sign-invariant CAD for {fi,f2,g) has 
14-87 cells using Qepcad 0/. We could then test a sample 
point of each cell to identify the ones where (j) is true. 

It is preferable to use the presence of ECs. Declaring an 
EC to Qepcad will ensure it uses the algorithm in f24f based 
on a single use of Pe{A) followed by P. Either choice re¬ 
sults in 289 cells. In particular, the solution set is described 
using 8 cells: all have y = 0 , z = —x but the x-coordinate 
unnecessarily splits cells at |(l±%/6). This is identified due 
to the projection polynomial d = disCy(res 2 (/i, g)). 

If we declare both ECs to Qepcad then it will use the 
algorithm in applying Pe{A) twice (allowed since n = 3) 
to produce a CAD with 133 cells. The solution set is now 
described using only 4 cells (the minimum possible). Note 
that d was no longer produced as a projection polynomial. 

For problems with n > 3 it is still possible to make use 
of multiple ECs. However, we must include the extra infor¬ 
mation necessary to provide order-invariance of the non-EC 
polynomials in the sections of ECs. The following theorem 
may be used to conclude that Pe{A) is appropriate. 

Theorem 3 ( [25] ). Let f and g be integral polynomials 

with mvar Xn, r = res(/, g), d = disc{g), and suppose r,d ^ 
0. Let S be a connected subset on which f is analytic 

delineable, g is not nullified and r and d are order-invariant. 

Then g is order-invariant in each section of f over S. 

Suppose we have a formula with two ECs, one with mvar 
Xn and the other with mvar Xn-i. The second could be ex¬ 
plicit in the formula or implicit (a resultant as described ear¬ 
lier). Theorem]^ allows us to use a reduced operator twice. 
We first calculate A' — Pe(A) where E contains the defin¬ 
ing polynomial of the first EC, and then Pf;,{A') where E' 
contains the defining polynomial of the other. Subsequent 
projections simply use P. When lifting we use Theorem 


to verify the CAD of R"“^ as order-invariant for Pf;i{A'); 
Theorem to verify the CAD of R"“^ order-invariant for 
E' everywhere and Theorem to verify it order-invariant 
for A' in the sections of E'- and Theorem and to verify 
the CAD of R" order-invariant for E and sign-invariant for 
A in those cells that are both sections of E and E'. 

2.2 Reductions in the lifting phase 

The main contribution of the present paper is to realise 
that the theorems above also allow for significant savings in 
the lifting phase of CAD. However, to implement these we 
must discard two embedded principles of CAD: 

1. That the projection polynomials are a fixed set. 

2. That the invariance structure of the final CAD can be 
expressed in terms of sign-invariance of polynomials. 

Abandoning the first is key to recent work in [11[ , while 

the second was also investigated in [10[ . 

2.2.1 Minimising the polynomials when lifting 

Consider Theorem it allows us to conclude that g is 
sign-invariant in the sections of / produced over a CAD 
of R"“^ order-invariant for 'P{/}({/.5})- Therefore, it is 
sufficient to perform the final lift with respect to / only 
(decompose cylinders according to the roots of / but not 
g). The decomposition imposes sign-invariance for / while 
Theoremguarantees it for g in the cells where it matters. 

Example 2. We return to Example^ Recall that desig¬ 
nating either EC and using f24f produced a CAD with 289 
cells. If we follow this approach but lift only with respect 
to the designated EC at the final step (implemented in our 
Maple package m) we obtain a CAD with I 4 I cells. 

This improved lifting follows from the theorems in [24| , 
but was only noticed 15 years later during the generalisa¬ 
tion of to the case of multiple formulae in [^ |^ . Exper¬ 
iments there demonstrated its importance, particularly for 
problems with many constraints (see Section 8.3 of [^). 

When we apply a reduced operator at two levels then we 
can make such reductions at both the corresponding lifts. 

Example 3. We return to the problem from Example\^ 
Set A = {/i, f 2 ,g} and E = {/i}. Then project out z using 

Pe{A) = {y^,y"‘ A 2xy^ + 2x^ A y^ - 1}. 

These are the resultants of fi with f 2 and g. The discrimi¬ 
nant of fi was a constant and so could be discarded, as was 
its leading coefficient (meaning no further coefficients were 
required). We set A' = Pe{A), E' = reSz(/i, / 2 ) = y^ and 

R — Tesy{y^, y^ -|- 2xy^ A 2x^ A y^ — 1) = {2x^ — 1)^. 

We have Pe'{A') = {i?} since the other possible entries (the 
discriminants and coefficients from E') are all constants. 
We hence build a 5 cell CAD of the real line with respect to 
the two real roots of R. We then lift above each cell with 
respect to y^ only, in each case splitting the cylinder into 
three cells about y — 0, to give a CAD o/R^ with 15 cells. 

Finally, we lift over each of these 15 cells with respect 
to fi to give 45 cells of R®. This compares to 133 from 
Qepcad, which used reduced projection but then lifted with 
all projection polynomials. No polynomials were nullified, so 
using Theorems^and^ the output is truth-invariant for (f. 






The additional lifting that Qepcad performed does not 
provide any further structure. For example, if we had lifted 
with respect to /2 at the hnal stage in Example then we 
would be doing so without the knowledge that it is deline- 
able. Hence splitting the cylinder at the sample point offers 
no guarantee that the cells produced are sign-invariant away 
from that point. So the extra work does not allow us to con¬ 
clude that /2 is sign-invariant (except on sections of /i). 

Note that using fewer projection polynomials for lifting 
not only decreases output size (and computation time) but 
also the risk of failure from non well-oriented input: we only 
need worry about nullification of polynomials we lift with. 

2.2.2 Minimising the cells for stack generation 

We can achieve still more savings from the theory in Sec¬ 
tion |2.1| by abandoning the aim of producing a CAD sign- 
invariant with respect to any polynomial, instead insisting 
only on truth-invariance for the formula. We may then lift 
trivially to cylinders over cells already known to be false, 
only identifying sections of projection polynomials if there 
is a possibility the formula may be true. The idea of avoiding 
computations over false cells was presented in [^. Our con¬ 
tribution is to explain how such cells can easily be identified 
in the presence of ECs. We demonstrate with our example. 

Example 4. Return to the problem from Examples^ — ^ 
and in particular the CAD of 9? produced with 15 cells in 
Example^ On 5 of these 15 cells the polynomial R is zero 
and on the others it is either positive or negative throughout. 

Now, 4> can only be satisfied above the 5 cells, as elsewhere 
the two EC defining polynomials cannot share a root and thus 
vanish together. We can already conclude the truth value for 
the 10 cells (false) and thus we do not need to lift over them, 
except in the trivial sense of extending them to a cylinder in 
R®. Lifting over the 5 cells where R = 0 with respect to fi 
gives 15 cells, which combined with the 10 cylinders gives a 
CAD o/R® with 25 cells that is truth-invariant for (f. 

This 25 cell CAD is not sign-invariant for f. The cylin¬ 
ders above the 10 cells in R^ where R f 0 may have fi 
varying sign, but since fi can never equal zero at the same 
time as /2 in these cells it does not affect the truth of (f. 

Identifying the 5 cells where i? = 0 in the CAD of R^ 
was trivial since they are simply the sections of the second 
lift, and hence those cells with second entry even in the cell 
index. Those sections produced in the third lift are similarly 
all cells where /i is zero, however, we cannot conclude that 
/2 is also zero on these. Theorem]^ only guarantees that /2 
is sign-invariant on such cells, so to determine those signs 
we must still evaluate the polynomials at the sample point. 

Reducing the number of cells for stack generation clearly 
decreases output size, and since the cells can be identihed us¬ 
ing only a parity check on an integer, computation time de¬ 
creases also. As with the improvements in Section r2.2.1[ this 
also decreases the risk of non well-oriented input: we only 
need worry about nullification over these identihed cells. 

3. ALGORITHM 

We present Algorithm to build a truth-invariant CAD 
for a formula in the presence of multiple ECs. We assume 
that the ECs are already identihed as input to the algorithm 
(they may have been hrst computed through propagation as 
described in Section . We assume further that each EC 
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Algorithm 1: CAD using multiple ECs 
Input : A formula <j) in variables xi,..., Xn, and a 

sequence of sets {Ek}^^!', each either empty or 
containing a single primitive polynomial with 
mvar Xk which dehnes an EC for 4>. 

Output: Either: D, a truth-invariant CAD of R" for (f 
(described by lists I and S of cell indices and 
sample points); or FAIL, if not well-oriented. 

Extract from (f the set of dehning polynomials 

for k = n,... ,2 do 

Set Bk to the hnest squarefree basis for prim(Afc); 
Set C to cont(Afc); 

Set Fk to the hnest squarefree basis for Ek; 
if Fk is empty then 

Set Afe_i := C'UP(Bfe); 
else 

it k = n or k = 2 then 
I Set Ak-i ■.= C U PFiiBi)- 
else 

L Set Ak-i ~CuPf{B,y, 

If El is not empty then set p to be its element; 
otherwise set p to the product of polynomials in Ai ; 
Build Di (Ii, ^i) according to the real roots of p; 

if n = 1 then 
^ return Di; 
for k — 2,... ,n do 

Initialise Dk = {lk,Sk) with Ik and Sk empty sets; 
if Fk is empty then 
Set L := Bk', 
else 

L Set L ~ Fk-, 
if Ek-i is empty then 

Set Ca := T>fc-i and Cb empty; 

else 

Set Ca to be cells in Dk-i with 1] even; 

Set Cb := Dk-i \ Ca', 

for each cell c £ Ca do 

if An element of L is nullified over c then 
return FAIL; 

Generate a stack over c with respect to the 
polynomials in L, adding cell indices and sample 
points to Ik and Sk', 
for each cell c £ Cb do 

Extend to a single cell in R*^ (cylinder over c), 
adding index and sample point to Ik and Sk', 

return = {I„,Sn). 


is primitive, and that all the ECs have diherent mvar (so in 
practice a choice of designation may have been made). 

Steps run the projection phase of the algorithm. 

Each projection starts by identifying contents and primitive 
parts. When there is no declared EC {Ei is empty) the pro¬ 
jection operator (T| is used (step[^. Otherwise the operator 
is used (ste p|l2[ ), unless it is the very hrst or very last 
projection (step |10p when we use f2f. In each case the out¬ 
put of the projection operator is combined with the contents 
to form the next layer of projection polynomials. 














Steps construct a CAD for the real line (and 

return it if the input was univariate). This is sometimes 
referred to in the literature as the base phase. If there is 
a declared EC in the smallest variable then the real line is 
decomposed according to its roots, otherwise according to 
the roots of all the univariate projection polynomials. 

Steps[^ —l^run the lifting phase, incrementally building 
CADs of R* for k = 2,..., n. For each k there are two 
considerations. First, whether there is a declared EC with 
mvar Xk- If so we lift only with respect to this (step 221 
and if not we use all projection polynomials with mvar Xk 
(step |20[ ). Second, whether there is a declared EC with mvar 
Xk-i- If so we restrict stack generation to those cells where 
the EC was satished. These are simply those with 1] 

(last entry in the cell index) even (step |26| |. We lift the other 
cells trivially to a cylinder in step |33| 

Algorithm clearly terminates. We will verify that it 
produces a truth-invariant CAD for the formula so long as 
the input is well-oriented, as defined below. 


Definition 2. For k = 2,... ,n define sets: 

• Lk — the lifting polynomials; the defining polynomial 
of the declared EC with mvar Xk if one exists, or all 
projection polynomials with mvar Xk otherwise. 

• Ck — the lifting cells; those cells in the CAD ofM.’‘~^ 
in which the designated EC with mvar Xk-i vanishes 
if it exists, and all cells in that CAD otherwise. 

The input of Algorithm^is well-oriented if for k = 2,...,n 
no element of Lk is nullified over an element of Ck. 


• When Ek is empty we use Theorem to assert all pro¬ 
jection polynomials with mvar Xk are order-invariant 
in the stacks over admissible cells giving (a) and (b). 

• When Ek is not empty and k = 2 we used the projec¬ 
tion operator §. Theorem allows us to conclude 
(b) and that the EC is sign-invariant in admissible 
cells. The stronger property of order-invariance fol¬ 
lows automatically since the lifting polynomials form 
a squarefree basis in two variables. 

• When Ek is not empty and fe = n we used the pro¬ 
jection operator @. Theorem [fallows us to conclude 
(b), but also (a) since in the case k — n the statement 
requires only sign-invariance. 

• When Ek is not empty and 2 < k < n we used the 
projection operator §. Theorem then allows us to 
conclude the statement. 

In each case the assumptions of the theorems are met by the 
inductive hypothesis, exactly over admissible cells as defined 
according to whether Ek-i was empty or not. 

From the dehnition of admissible cells, we know that (f> 
is false (and thus trivially truth invariant) upon all cells 
in the CAD of R" built over an inadmissible cell of R*’, 
k < n. Coupled with the truth of (a) for k — 1,..., n, this 
implies the CAD of R" is truth-invariant for the conjunction 
of ECs (although it may not be truth-invariant for any one 
individually). The truth of (b) implies that on those cells 
where all ECs are satisfied, the other polynomials in cf are 
sign-invariant and thus ij) is truth-invariant. □ 


Theorem 4. Algorithm^ satisfies its specification. 

Proof. We must show the CAD is truth-invariant for 
unless the input is not well-oriented when FAIL is returned. 

First consider the case where n — 1. The projection phase 
would not run, with the algorithm jumping to the CAD 
construction in step |13[ returning the output in step |16| If 
there was no declared EC then the CAD is sign-invariant 
for all polynomials defining cp and thus every cell is truth 
invariant for (p. If there was a declared EC then the output 
is sign-invariant for its defining polynomial. Cells would 
either be intervals where the formula must be false; or points, 
where the EC is satished, and the formula either identically 
true or false depending on the signs of the other polynomials. 

Next suppose that the input were not well-oriented (Def¬ 
inition]^. For a hxed k, the conditional in steps [T^ — |22| 
sets the lifting polynomials Lk to L and the conditional in 
steps — [^the lifting cells Ck to Ca- Thus it is exactly 
the conditions of Dehnition]^ which are checked by step |29[ 
returning FAIL in step[^when they are not satished. If the 
lifting phase completes then the input is well-oriented. 

From now on we suppose n > 1 and the input is well- 
oriented. For a hxed k dehne admissible cells to be those 
in the induced CAD of R*’”^ where all declared ECs with 
mvar smaller than Xk are satished, or to be all cells in that 
induced CAD if there are no such ECs. Then let I (fe) be the 
following statement for the CADs produced by Algorithm]^ 
Over admissible cells (in the CAD o/R*^ is: 

(a) order-invariant for any EC with mvar Xk; 

(b) order- (sign- if k = n) invariant for all projection poly¬ 
nomials with mvar Xk on sections of the EC over admis¬ 
sible cells, or over all admissible cells if no EC exists. 

We have already proved /(I), and I{n) may be proved by 
induction. To assert the truth of /(fe) we note the following: 


4. WORKED EXAMPLE 

Assume variable ordering z)~y>~x)~u>~v and dehne 
/i ■.= X — y z^, f2 := -b — 1, g := x^ — 1, 

/a := a;-b y-b fi-.= z'^-\-u^-v'^-l, h\=z. 

We consider the formula 

<() = /i=0A/2 = 0A/3 = 0A/4 = 0A(?>0A/i>0. 

The solution can be found manually by decomposing the 
system into blocks. The surfaces /i and /a are graphed in 
(a;, y, z)-space on the left of Figure They meet only on 
the plane y = 0 and this projection is shown on the right. 
The surfaces /a and /4 are graphed in {z, u, u)-space on the 
left of Figure]^ and meet only when 2 = ±1. We consider 
only 2 = -bl due to fe > 0, with this projection plotted on 
the right. We thus see that the solution set is given by 

{u = ±u, X = —1, y = 0, 2 = 1}. 

To ascertain this by Algorithm we must hrst propagate 
and designate ECs. We choose to use /i hrst, calculate 

res 2 (/i,/ 2 ) = {-u +v'^ -x + y-if 

and assign ri to be the square root: the dehning polynomial 
for an EC with mvar y. Similarly consider 

resy (ri,resz(/i, /a)) = 16(u^ - -b a; -b 1)"^, 

resy (ri,reSz(/i,/4)) = A{u - v'^ f 

and assign ra := -b a; -b 1, ra := — v^: dehning 

polynomials for ECs with mvar x and u respectively. There 
is no series of resultants that leads to an EC with mvar 
u. We hence have := {/i}, {ri}, {ra}, {ra}, { } as 

input for Algorithm]^ along with (p. 



Figure 2: The polynomials fi and /a from Section]^ 


Figure 3: The polynomials /a and /4 from Section]^ 


The algorithm starts by extracting the defining polynomi¬ 
als A5 = {/i, / 2 , /a, fi, g, h} and finds B5 = A5,F5 = E5 (in 
fact Fi = Ei for alH = 1,..., 5). There is a declared EC for 
the first projection so we use the operator § to derive 

A 4 ;= Pp^iBs) = {(x^ - if, {-u^ - X Ay - if, 

{u - - X A y - if ,x - y}. 

Hence C := {x^ — 1} and 

Bi ~ {y, y - X, -u^ A - x A y - l,u^ - - x A y - 1]■ 

For the next projection we must use operator § , giving 

A-i := CVJPpfBi) = {x^ — l,u^—v^AxAl,u^—v^,u^—v^Al} 

noting that for this example the extra discriminants in 
all evaluated to constants and so could be discarded. Then 

B3 := {a;^ ——C := ,u^ — Al}, 

and the next projection also uses to produce 

A _ f 2 2 2 2 , -1 4 r>22| 4|r>2 riSi 

A 2 := {u — V ,u — V -I- 1,u — 2u V A v A 2u — 2v 

For the final projection there is no EC and so we use operator 
Q to find Ai ;= {u^}. The base phase of the algorithm 
hence produces a 3-cell CAD of the real line isolating 0. 

For the first lift we have L = {u^ ~ and Ca containing 
all 3 cells. Above the two intervals we split into 5 cells by the 
curves u = ±u, while above u = 0 we split into three cells 
about the origin. From these 13 cells of we select the 5 
which were sections of for Ca- These are lifted with 

respect to L = {ra}, and the other 8 are simply extended to 
cylinders in Together this gives a CAD of R^ with 23 
cells. The next two lifts are similar, producing first a CAD 
of R"^ with 53 cells and finally a CAD of R® with 113 cells. 
The entire calculation takes less than a second in Maple. 

Choice in EC designation 

Algorithm could have been initialised with alternative EC 
designations. There were the 4 explicit ECs with mvar z. 


and by taking repeated resultants we discover the following 
implicit ECs, organised in sets with decreasing mvar: 

r22 2| I't 2|2| i'i 

{y ,u —V Ax — yAl,—u -I-u A x — y A I, 

A X A y A I, —u^ Av^AxAyAl}, 

{x A 1, —u^ -I- -I- 2: -I- 1, -I- 2 ; -I- 1}, {u^ — v^}. 

There are hence 60 possible permutations of EC designation, 
but they lead to only 3 different outputs, with 113, 103 and 
93 cells. Heuristics for other questions of CAD problem 
formulation [15[ |21[ |31| could likely be adapted to assist 

here. We note that 93 cells is not a minimal truth invariant 
CAD for as it splits the CAD of R^ at u = 0 (identified 
from the discriminant of the only EC with mvar u). 

Comparison with other CAD implementations 

A sign-invariant CAD of R® for the 6 polynomials in the 
example could be produced by Qepcad with 1,118,205 cells. 
Neither the RegularChains Library in Maple nor our 
Maple package could produce one in under an hour. 

Our implementation of [^, which uses operator § once 
but also performs the final lift with respect to the EC only, 
can produce a CAD with either 3023, 10935 or 48299 (twice) 
cells depending on which EC is designated. The Qepcad 
implementation of [24| gives 11961, 30233, 158475, or 158451 
cells. Comparing these sets of figures we see the dramatic 
improvements from just a single reduced lift. 

Allowing Qepcad to propagate the 4 ECs (so a similar 
projection phase as Algorithm but then a normal CAD 
lifting phase) produces a CAD with 21079 cells. By declar¬ 
ing only a subset of the 4 (which presumably changes the 
designations of implicit ECs) a CAD with 5633 cells can be 
produced, still much more than using Algorithm 

The RegularChains Library can also make use of multiple 
ECs, as detailed in [^. The version in Maple 18 times out 
after an hour, however, with the development version a CAD 
can be produced instantly. There are choices (with analogies 
to designation [^) but they all lead to a 137 cell output. 
In particular, they all have an induced CAD of the real line 
which splits at u = ±1 as well as u = 0. 

5. COMPLEXITY ANALYSIS 

We build on recent work in to measure the dominant 
term in bounds on the number of CAD cells produced. Nu¬ 
merous studies have shown this to be closely correlated to 
the computation time. We assume input with m polynomi¬ 
als of maximum degree d in any one of n variables. 

Definition 3. Consider a set of polynomials pj. The 
combined degree of the set is the maximum degree (taken 
with respect to each variable) of the product of all the poly¬ 
nomials in the set: maxi(deg,^. ( Pj ))• 

The set has the (m,d)-property if it may be partitioned 
into m subsets, each with maximum combined degree d. 

For example, {y^ — x,y^ A 1} has combined degree 4 and 
thus the (1, 4)-property, but also the (2, 2)-property. 

This property (introduced in McCallum’s thesis) can mea¬ 
sure growth in the projection phase. In we proved that 
if A has the (m, d)-property then P{A) U cont(A) has the 
(M, 2d^)-property with M = (|(m-|- 1)^J. When m > 1, 
we can bound M by (but we need 2m^ to cover m = 1). 











If A has the (m, d)-property then so does its squarefree 
basis. Hence applying this result recursively (as in Table 
measures the growth in (m, ci)-property during projection 
under operator 0 - After the first projection there are mul¬ 
tiple polynomials and so the tighter bound for M is used. 

The number of real roots in a set with the (m, dj-property 
is at most md. The number of cells in the CAD of is thus 
bounded by twice the product of the hnal two entries, plus 
1. Similarly, the total number of cells in the CAD of R" by 

{2md + 1) nrii [2 + l] . (4) 

Omitting the -I-Is will leave us with the dominant term of 
the bound, which evaluates to give the following result. 

Theorem 5. The dominant term in the bound on the 
number of CAD cells in R" produced using 0 is 

. (5) 

From now on assume I ECs, 0 < ^ < min(m,n), all with 
different mvar. For simplicity we assume these variables are 
x„-£+i (the first £ projections are reduced). 

Lemma 6. Suppose A is a set with the {m,d)-property 
and E C A has the {l,d)-property. Then cont(A) U Pf;{A) 
has the {2m, 2d^)-property. 

Proof. In we proved that applying Pe{A) U cont(A) 
gives a set of [^(3m + 1 )J polynomials of combined degree 
2df. The extra m — 1 discriminants required by operator 0 
will each have degree at most d{d — 1 ), so pairing them we 
have [|(jn — 1)] sets of combined degree at most 2d?. Then 

L|(3m-bl)J -brl(m-l)] =m + L|(m-bl)J -f [f J 

and since m £ Z this always equals 2m. □ 

We apply this recursively in the top half of Table with 
the bottom derived via the process for P, as in Table 

Define di and m; as the entries in the Number and Degree 
columns of Table |2] from the row with i Variables. We can 
bound the number of real roots of projection polynomials 
in i variables by midi. If we lifted with respect to all these 
projection polynomials, the cell count would be bounded by 

n"=r [2m,d, + 1] = ULo [2 {rm2^‘-^d^^) + l] 

• [2 + 1] . (6) 

Omitting the -fl from each product allows us to calculate 
the dominant term of the bound explicitly as 

(7) 

Now we consider the benefit of improved lifting. Start by 
considering the CAD of There can be no reduced 

lifting until this point and so the cell count bound is given 
by the second product in 0 , which we will denote by j. 
The lift to R"“^ will involve stack generation over all cells, 
but only with respect to the EC. This can have at most d„-t 
real roots and so the CAD at most [2dn-e -I- l](t) cells. 

The next lift, to will lift the sections with respect 

to the EC, and the sectors only trivially (to produce the 
same number of cylinders). Hence the cell count bound is 
[2dn-(e-i) -\- l]dn-i{]) -\- {dn-i -\- l)(t) with dominant term 
2dn-(e-i)dn-e{'\)- Subsequent lifts follow the same pattern 
and so 2 d„d„_i .. . d„_(f_i)d„_f (|) is the dominant term in 
the bound for R". This evaluates to give the following result. 


Table 1: Projection under operator (Q. 


Variables 

Number 

Degree 

n 

m 

d 

n — 1 

2m^ 

2d^ 

n — 2 

4m^ 

• 

00 

n — r 



1 

r\TL — 2 oTT.— 1 

2^ 


Table 2: Projection with ([3|) 

Variables Number 

£ times and then 

Degree 

n 

m 

d 

n — 1 

2m 

2d^ 

n — £ 

n — {£-\-l) 

2^m 

2^^m? 

2^'-^d?'' 

22^+1-1^2^+1 

n — {£-\-r) 


22<^+--1^2'^+'- 

1 2^ 

^n- 



Theorem 7. Consider the CAD of M." produced using Al- 
gorithm^in the presence of ECs in the top £ variables. The 
dominant term in the bound on the number of cells is 


2nLo [2'*"“^d2“] [2 (2'^"^m^"2'^'*"-^d^ 


= (2a) m 2 


)] 

( 8 ) 


The bound in Theorem is strictly less than the one in 
Theoremj^ The double exponent of m has decreased by the 
number of ECs; the result of the improved projection in 0 - 
Improved lifting reduced the single exponents further still. 


6. CONCLUSIONS AND FUTURE WORK 

We have explained how the existing theory for CAD pro¬ 
jection using ECs can also be leveraged for significant sav¬ 
ings in the lifting phase. We can reduce both the projection 
polynomials used for lifting and the cells over which stacks 
are generated. We have formalised these ideas in Algorithm 
verified their use in Theorem and demonstrated the 
benefit with a worked example and complexity analysis. 

A key question is how to best deal with non-primitive 
ECs? Consider cp := zy — 0 A g}. Under ordering • • • 
z y y y ... the EC zy = 0 is not primitive, so Algorithm 
[^cannot use it. We may be tempted to take E = { 2 } as 
the primitive part, project with operator S and include 
the content y in the first projection. The CAD of {y,... )- 
space would be sign-invariant for y and thus the CAD of 
{z, y,... )-space truth invariant for the EC (over admissible 
cells). But we can no longer say only sections are admissible 
for the next lift as there may be cells with 2 7 ^ 0 and y = Q. 
We could instead lift over all cells. Alternatively we might 
rewrite (p as (p := {z = 0A<p)V {y = 0A(^), so each clause has 
its own EC. The theory of truth-table invariant CADs 
is designed to deal with such input, but would require its 
own extension to use beyond the first projection. Of course, 
this extension would also be valuable in its own right. 
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